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I. INTRODUCTION 

Metallic single-electron transistors (SETs) have found 
several applications including charge detectors and am- 
plifiersii^, electron pumps for metrology^ii^ and low- 
temperature thermometers-. More recently, single 
dopants in silicon have emerged as promising candidates 
for qubits of the ultimately small size^. Silicon tech- 
nology offers the opportunity both to realize classical 
SETsi"— and to probe single dopantsi^""— . The use of 
the SET for quantum-mechanical measurement of the 
dopant qubit requires characterization of the coupling 
strength between SET and dopant atom. Here we study 
the mutual capacitive and tunnel couplings between a 
dopant atom and a SET, realized in a silicon nanowire 
field-effect transistor. 

The interplay between Coulomb interaction and 
tunneling is most thoroughly studied in mesoscopic 
physics^"', where it leads to the so-called Coulomb block- 
ade (CB) oscillationsSiS of conductance, commonly ob- 
served in metallic SETs^^ and semiconductor quantum 
dots— The energy scale governing the CB oscilla- 
tions in a metallic SET is set by the charging energy 
Ec ~ e^/2C, where C is the total capacitance of the 
SET island relative to the surrounding conductors, and 
— e is the electron charge. At temperatures T ~ Ec, the 
conductance as a function of the gate voltage acquires 
an oscillatory component and, with lowering the tem- 
perature, it turns into a series of spikes (CB peaks), at 
T ^ Ec- At the lowest temperatures, the height of a 
CB peak depends on the properties of the electron wave 
function inside the SET; namely, on its values at the con- 
tacts. The peak height changes randomly from peak to 
peak, due to the chaotic nature of the electron motion 
(and wave functions) inside the SET. Such mesoscopic 
fluctuations of the CB peaks persist up to temperatures 
on the order of the single-particle level spacing, AE, of 
the SET. It is convenient to distinguish between SETs 
and quantum dots by the magnitude of AE: SETs are 
devices in which AE <^ Ec and an intermediate tem- 
perature regime, AE T ^ Ec, is possible, whereas 
quantum dots are devices in which AE ^ Ec- The exis- 



tence of the temperature regime AE <C T <C Ec makes 
SETs outstanding candidates for applications. 

In the temperature regime AE ^ T ^ Ec, the meso- 
scopic fluctuations of the CB peaks are suppressed. The 
CB peaks have equal height, which is determined by the 
transparencies of the SET tunnel barriers, and thus, it is 
independent on the details of the electron motion inside 
the SET island. The peak height remains constant over 
a large number of CB oscillations and may vary weakly 
with the gate voltage due to the way the gate couples 
to the SET tunnel barriers. A strong modulation of the 
peak height on a scale of several CB oscillations is usually 
considered as an anomaly in this regime. Such anoma- 
lies in the CB oscillations have been attributed to charge 
traps, which may be present in the SET surrounding2^i2£. 
Charge traps is a common and, typically, unwanted in- 
gredient in modern electronic nanodevices. Their most 
conspicuous collective effect is the voltage switch, which 
makes devices difficult to control. In silicon nanostruc- 
tures based on doping modulation, charge traps occur 
naturally at the borders of the doped regions, where the 
dopants are likely to diffuse away from their majority 
distribution and form nearly isolated charge traps. In 
this paper, we study the effects associated with a single 
charge trap (a dopant atom) in the CB oscillations of a 
SET. 

We observe an interplay between mutual capacitive 
coupling and tunneling in the silicon SET coupled to a 
dopant atom. In particular, we observe a spectacular en- 
hancement of the conductance through the SET when 
transport occurs by resonant tunneling via the dopant 
atom. We develop a theory which takes into account both 
tunneling and capacitive coupling between dopant atom 
and SET, and assess the coupling strengths by analyzing 
the CB oscillations of the linear conductance. The mu- 
tual capacitive coupling between SET and dopant atom, 
extracted from the experiment, has an energy scale of 
^ 1 meV, which indicates a good charge sensitivity of 
the SET. In this double-dot system, i.e. the system of 
coupled dopant atom and SET, the SET can be used for 
manipulation and readout of the single dopant atom. 

The paper is organized as follows. In Sec.|lTl we present 
the experimental results alongside with a qualitative ex- 
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planation of the observed effects and a quantitative fit 
to the experimental data using the theory developed in 
subsequent sections. We also formulate a simple elec- 
trostatic model which helps in discussing the experimen- 
tal findings. In Sec. IIIIl we proceed with a systematic 
analysis of the electrostatic model. We first consider the 
low-temperature limit, see Sec. IIII Al in which transport 
is activational, and hence, governed by the energy spec- 
trum derived from the electrostatic model. Then, we 
complete the electrostatic model by kinetic-energy and 
tunneling terms and proceed with the study transport 
through the coupled donor-SET system by means of rate 
equations, see Sec. IIIIBI We find that at low tempera- 
tures, the linear conductance could also be derived from 
a circuit approach, in which the donor-SET system is re- 
placed by a circuit of resistors. We analyze the validity of 
this approximation at higher temperature by solving the 
rate equations numerically and find good agreement for 
the parameter sets which we extracted from fitting to the 
experimental data. In Appendix|Xl we give the necessary 
formulae for a quantitative description of the double-dot 
stability diagram occurring in the electrostatic model. 



II. MAIN RESULTS 

We implement an electrostatically-defined SET in a 
Si nanowire fabricated using the industrial etching and 
doping techniques described in Ref . [l^- A schematic view 
of the device is shown in Fig. [T^. The SET island is 
formed by electron accumulation under a top gate in a 
region of intrinsic Si, see Fig. [T|d. The leads are defined 
by doping with As atoms (donors). 40- nm- wide spacers, 
made of silicon nitride, are deposited on both sides of 
the gate prior to the ion implantation in order to create 
low-doped regions acting as tunnel barriers between the 
leads and the island, see Figs.[T^ and[T]D. Indeed, samples 
without spacers do not exhibit SET behavior, but allow 
to study transport through single dopants which diffused 
away from the leads into the Si region below the gateii. 
A large number of devices with spacers 10^) have 
been fabricated and more than 50 have been studied at 
low temperatures, T > 40 mK. All devices exhibit low 
1/f noise and very periodic CB oscillations. 

The high quality of the SET, which can be inferred 
from Figs. [T]: and [TJl, allows us to focus on special re- 
gions of gate voltage, in the trace of CB osciUations, 
where the transport properties of the SET are strongly 
perturbed by single As atoms. We observe between 1 and 
5 of such anomalies in each device in different ranges of 
the top-gate voltage Vg. In previous reports of anoma- 
lous behaviour of the SET— i^i^— , the coupling between 
the SET island and the single dopants has been predom- 
inantly capacitive. Here we observe dramatic increase of 
the current through the SET due to resonant tunneling 
through single dopants. In Fig. [21 we show the signatures 
of tunnel coupling of an As atom to the lead and SET 
island; resonant features are observed both in the linear 



and differential conductances. 

We focus on the linear conductance and analyze two 
representative types of anomalous behaviour observed in 
the experiment, see Figs. [3^ and[5}D. The anomalous be- 
haviour has the following generic pattern: (i) at the low- 
est temperatures, the CB peaks are strongly suppressed 
within an interval of Vg, (ii) at intermediate tempera- 
tures, two maxima develop in the envelope of the CB os- 
cillations on both sides of the interval of suppression, (iii) 
at higher temperatures, the region of suppressed conduc- 
tance turns into a region of elevated conductance, with 
a single maximum in the envelope. This pattern, evi- 
denced by Fig. [31d, is generic to the case when both ca- 
pacitive and tunnel couplings are present. The measure- 
ments shown in Fig. are consistent with having only 
capacitive coupling. In this case, the envelope of CB os- 
cillations exhibits no maxima; the anomalous behaviour 
consists in (i) only. 

We successfully reproduce the anomalous behaviour 
of the CB oscillations observed in the experiment, tak- 
ing into account capacitive and tunnel couplings of the 
dopant atom to the SET. The results of our theory are 
shown in Figs. |3t and |3Ji, and correspond, respectively, 
to Figs. [3^ and |3)d of the experiment. In the following 
paragraphs, we formulate a simple model for the elec- 
trostatic energy of the device and explain the observed 
behavior in terms of sequential tunneling of electrons be- 
tween energetically favorable states. 

The dopant atom (donor) resides in one of the SET 
barriers; let it be the barrier next to the source. A sketch 
of the setup is shown in Fig. 2^. The energy associ- 
ated with the charges in the device depends on the donor 
charge n = 0, 1 and the SET charge N (both measured in 
units of the electron charge, — e), and can be abbreviated 
in the following form 

E{n, N) = ean + Ec{N - Ngf + Ui2n{N - Ng), (1) 

where is the energy of the donor level, Ec is the charg- 
ing energy of the SET, Ng — CgVg/e is the dimensionless 
gate voltage, with Cg being the gate-to-SET capacitance, 
and Ui2 is the energy of mutual capacitive coupling. The 
top gate couples to both the SET and the donor. How- 
ever, the lever arm of the donor is considerably smaller 
than that of the SET, because the donor is more dis- 
tant from the gate than the SET island is, and the donor 
charge is partly screened due to proximity to the lead. We 
express this fact in the relation = —aEcNg + const, 
with a 1. 

At low temperatures, transport occurs due to tran- 
sitions from the ground state to the lowest-in-energy 
charge configurations. The stability diagram, commonly 
used for double dots'^'^, shows the ground state as a func- 
tion of and Ng, see Fig.|l|D and Appendix El Variation 
of the gate voltage Vg leads to traversing the plane of the 
stability diagram along an inclined line — the "work- 
ing line" of the device. In Fig. Hb, the dash-dotted line 
represents the working line. Each time the working line 
crosses a solid line in the stability diagram, two charge 
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FIG. 1. Gated nanowire device, (a) Geometrical arrangement of the main components of the device. The Si nanowire channel 
(20 nm thick and 60 nm wide) bridges the constriction between the source and drain. The gate (40-nm-wide poly-silicon finger) 
is deposited on top of the nanowire and it is surrounded by silicon nitride (Si3N4), thus forming two spacers (each 40 nm wide) 
on both sides of the gate. The red spheres illustrate the doping with As atoms inside the Si material; dopants on the surface 
of Si are omitted for clarity, (b) Cross section along the nanowire. In the experiment, the implantation with As atoms is made 
after the deposition of gate and spacers. Therefore, during the implantation, a small number of As atoms may diffuse into the 
region under the spacer, forming isolated donors. The SET island is formed by electron accumulation when a positive voltage 
is applied to the gate, (c) Transmission-electron-microscopy image of a cross section along the nanowire in an actual device. 
The nanowire (Si) is separated from the poly-silicon finger (gate) by a 4nm thick oxide layer (Si02); the spacers are hard to 
discern on the micrograph, (d) High-resolution image showing the atomic lattice of the nanowire and its interface with the 
buried oxide (BOX). 
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FIG. 2. Experimental signatures of tunnel coupling between a single dopant atom and a single-electron transistor, (a) Linear 
conductance through the silicon nanowire as a function of the gate voltage Vg at various temperatures. The strong, resonance- 
like modulation of the CB oscillations is attributed to the presence of a dopant atom in the SET barrier, (b) Color plot of the 
differential conductance of the device versus the gate and bias voltages. The small Coulomb diamonds of the SET are strongly 
modulated in intensity, revealing, on a larger scale, resonant tunneling associated with the dopant atom. 



configurations have equal energies in Eq. ([T|). (For sim- from the Fermi level in the leads.) The conductance, 
plicity of argument, we measure single-particle energies however, may display CB peaks only for certain types of 
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degeneracies. 

We consider first the effect of direct tunneling tfirougli 
the SET, which does not involve a change of the donor 
charge. The activation energy of this process is deter- 
mined by the difference E{n,N) — E{n,N + 1). Note 
that both terms in the latter expression are evaluated 
at the same n, which corresponds to the condition that 
there is no current through the donor. That energy dif- 
ference reaches zero every time the working line of the 
device intersects a vertical line of the stability diagram, 
which is possible for Ng < N~ and Ng > , see Fig.|3)D. 
These intersections correspond to CB peaks which have 
temperature-independent heights (zero activation energy, 
see Fig. lit). For Ng < Ng' , the conductance around the 
peaks can be approximated by^ 



G(7V5)~^g(x-a;Ar), g{z) ^ 



9l9r 



N 



h gL + 9r sinh(z) 



(2) 

with X = 2EcNg/T and xat = 2£'c(iV + l/2)/T. Here, 
the summation goes over integer N, the function g{z) 
describes the shape of a single peak, and gn are the 
dimensionless conductances, in units of 2e^//i, of the SET 
tunnel barriersSS,. Equation ([5]) describes transport due 
to sequential tunneling at low temperatures (T ^ Ec) 
and for positions in the stability diagram far above the 
triple points, where the donor degree of freedom is frozen 
to n = 0. For Ng > N^, far below the triple points, 
where the donor degree of freedom is frozen to n = 1, 
the conductance is given by Eq. ©, with xjv xn + 
U12/T. The periodic CB oscillations on both sides of 
the anomaly are shifted with respect to each other by 
a fraction, Ux^/^Ec, of the period. Using this fact, we 
extract C/12 = 1.2±0.1 meV for the sample in Fig.[3K and 
U12 — 1.15 ±0.1 meV for the sample in Fig. 

Charge-charge correlations are essential to understand- 
ing the transport mechanism in the anomaly region, 
Ng < Ng < N^ . The donor charge adjusts to the 
SET charge, lowering the electrostatic energy in Eq. ([1]). 
This effect is strongly pronounced around the SET res- 
onances, at which the difference E{n,N) — E{n,N + 1) 
vanishes, see dashed lines in inset of Fig. HJd. For con- 
creteness, let us discuss the dashed lines of the upper 
triple points, where E{0,N) = E{0,N + 1). The lowest- 
in-energy charge configuration at these points is (1,7V), 
whereas the configurations (0,7V) and (0,7V -I- 1) corre- 
spond to excited states. In order for the transport to 
occur, the donor charge has to change from n = 1 to 
n = 0. This process is suppressed by an activation en- 
ergy 6 = E{0,N) - E{1,N), evaluated at the SET res- 
onance E{0, 7V) — E{0, N + 1). The presence of a finite 
activation energy is responsible for the behavior of type 
(i) observed in the experiment. The height of a CB peak 
in the anomaly region decreases with lowering the tem- 
perature as oc exp(— (5/T), where S is shown in Fig. |4]; 
for the whole range of the gate voltage. Note that the 
gap 5 reaches maximum, S^^^x = (t^i2/2)(l — Ui2/2Ec), 
in the middle of the anomaly and it goes to zero towards 



the edges. We analyze the conductance around a CB 
peak in Sec. IIIII At low temperatures (T <C Ec), the 
conductance can be written as 



N 



-{n)Q{x - xjv - U12/T)] 



(3) 



where {n) is the average occupation of the donor. The oc- 
cupation (n) varies strongly with the gate voltage inside 
the anomaly region, whereas it approaches a constant 
value of or 1 outside. We give explicit expressions for 
(n) in the vicinity of triple points in Sec. IIII A3I Equa- 
tion ([3]) has a simple physical meaning: the donor charge 
shifts the gate voltage "seen" by the electron passing 
through the SET; the probabilities of having the donor 
charge 1 and are, respectively, (n) and 1 — (n). 

The CB peaks in the anomaly region described by 
Eq. ^ occur at shifted positions, TVg = 7V + 1/2 + 
where (/s is a phase shift of CB oscillations. We show the 
behaviour of (p at low and high temperatures in Fig. |4}i. 
At T C/12, the phase shift is given hy ip — {n)Ui2/2Ec, 
where (n) assumes the form of the average occupation of 
an isolated donor. In addition to the CB peaks being 
shifted in the anomaly region, they also become broader 
than usual. Usually, the width of a CB peak is given 
by ANg oc T/Ec, which is a relation commonly used 
in determining the electron temperature in the sample. 
We find that, a,t T < S, this relation is no longer valid 
and one has ANg oc nim{T/Ui2,d/Ec')- With the peaks 
being broader in the anomaly region, the conductance 
in the middle of CB valleys becomes larger. Indeed, in- 
stead of the usual suppression of the conductance in the 
middle of the valley, oc exp{—Ec/T), we find a smaller 
suppression factor, oc exp(— i5/T), with the activation en- 
ergy 6 < Ec reaching its minimum, Jmin = Ec — U12/2, 
in the middle of the anomaly. 

We concentrate now on the role of tunneling through 
the donor and, for that purpose, we dispense with the 
possibility of direct tunneling between the SET and the 
source {gL = 0). To have an activationless transport 
in that case one has to satisfy simultaneously two con- 
ditions for the energy: E{ti,N) — E{n,N + 1) and 
E{0, 7V+ 1) = E{1, N). Clearly these two conditions may 
be satisfied only at the triple points. We consider a small 
a for which the working line passes by triple points suffi- 
ciently close (a <C T/ Ec)- The maxima of the CB peaks 
form a smooth envelope function, Genv{Ng), which con- 
tains information about the tunnel coupling of the donor. 
Indeed, we find that, in the vicinity of the upper row of 
triple points. 



G.nANg) 
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h TiTL + Tn) cosh2(y/2)' 



(4) 



with y = aEc{Ng - Ng)/T and Tl {Tr) being the 
tunneling rate between donor and source (SET island). 
Equation Q coincides with the conductance through a 
resonant level. In deriving Eq. we made several 
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FIG. 3. Anomalies in transport through SETs caused by single dopant atoms, (a) Linear conductance measured in a sample 
in which the dopant atom couples mostly capacitively to the SET. At T = 1 K, the CB peaks are suppressed in an interval 
of Vg, marking the extension of the anomaly. At T = 4.2 K, the CB oscillations are weakly modulated, (b) Modulation of 
CB oscillations in a different device, in which the dopant atom is tunnel-coupled to both the SET island and the lead. The 
envelope of CB oscillations reveals resonant tunneling through the dopant atom. At T = 0.9 K, two resonances, separated by a 
suppression region, are visible in the envelope function. At T = 4.2 K, the resonances overlap and form a single broad resonance 
due to tunneling via the dopant atom. Note the factor-of-5 difference in the scale on the ordinate axis as compared to (a) . (c) 
Linear conductance calculated using a circuit approach, see Sec. IIII B 41 With the donor coupling parameters U12 ~ 1.2 meV 
and Tl = Fij = 0, the experiment in (a) is reproduced almost quantitatively, (d) The same as in (c), but plotted using the 
coupling parameters I7i2 = 1.15 meV and Vl ~ Lu — 40 /xeV, which nears the result to (b). The SET parameters used in (c) 
and (d) are tabulated in Table IH 



simplifying assumptions (see Sec. IIII ]T3| . which do not 
change the essence of our result: the information about 
the donor is encoded in the envelope of the CB oscilla- 
tions. At low temperatures, T ^ U12, the envelope func- 
tion consists of two resonances centered at Ng = N^, 
see Fig. The resonance at Ng = differs slightly 
from the one at Ng = Ng , because of the spin degen- 
eracy on the donor. Each of the resonances corresponds 
to traversing a horizontal line of triple points in Fig. 2)3. 
The behavior of type (ii) observed in the experiment (see 
Fig. inb) is, thus, due to resonant tunneling through the 
donor level. At high temperatures, T ^ U12, the two 
resonances overlap and merge into a single broad fea- 
ture, see Fig. HJj, corresponding to the behavior of type 
(iii) seen, e.g., in Fig. |3}3. 

Thus far, we have separated the effect of tunneling 
through the donor from direct tunneling by setting ql = 



0. In this special case, the CB peaks are suppressed at 
low temperatures both inside and outside the anomaly 
region. The corresponding activation energy d is shown 
by the dashed line in Fig. At a finite but small ql, 
the envelope maxima around Ng = N^ remain present, 
but the smaller activationless CB peaks are visible out- 
side the interval N~ < Ng < iV+ . The activation energy 
6 is, then, given by the solid line in Fig. In order to 
describe the interplay between direct tunneling and tun- 
neling via the donor, we employ a circuit approaches,. We 
associate a resistor with each tunnel junction and sum up 
the resistances of the so-obtained circuit according to the 
classical laws for circuits, see Sec. IIII B 41 The end result 
is a conductance formula, which we use to reproduce the 
experimental data in Fig. |31 We also back up the circuit 
approach by an exact numerical evaluation of the conduc- 
tance. We find that the circuit approach works extremely 
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FIG. 4. Setup, stability diagram, energy gap, phase shift, 
and envelope of CB oscillations, (a) Schematic view on the 
device. The SET island, with charge A^, is tunnel-coupled 
to two leads (source and drain). The donor, with charge n, 
is situated in the tunnel barrier between the source and the 
SET island. The top gate, at voltage Vg, couples capacitively 
to the SET island and donor; the coupling to the latter is the 
weakest of the two. (b) Stability diagram of the donor-SET 
system. Solid lines separate regions with different ground 
state charge configurations. Sweeping the gate voltage Vg in 
the experiment corresponds to traversing the plane of the sta- 
bility diagram along the dash-dotted line - the "working line" 
of the device. The anomaly occurs in the region between the 
upper and lower rows of triple points. For small coupling of 
the donor to the gate (a <C 1), the working line is nearly hor- 
izontal; the anomaly region comprises, then, a large number 
of charge configurations, (c) The energy gap, governing the 
activated conductance at CB peaks, as a function of the gate 
voltage Ng. The solid line corresponds to the case ql 0, 
whereas the dashed line to ql = 0. (d) The phase shift ip, 
which gives the position of the CB peaks across the anomaly 
region and outside. At low temperatures, ip changes stepwise 
in the center of the anomaly. At high temperatures, (p changes 
gradually over an interval which grows proportionally to the 
temperature, (e) The envelope of the CB oscillations versus 
the gate voltage. At low temperatures, two resonances oc- 
cur on the sides of the anomaly. At high temperatures, the 
resonances overlap, forming a single broad resonance. 



well for the experimentally relevant parameters. 



III. DERIVATION OF MAIN RESULTS 



presented in Sec. [IT] regarding our theoretical description 
of transport. 



Activation energies derived from the 
electrostatic model 



1. Activation energy at the CB peak 

With lowering the temperature, the SET conductance 
due to sequential tunneling is suppressed as 



Gcxexp(-A/T)^ 



(5) 



where A is an activation energy and T is the temper- 
ature. The activation energy is the excitation energy 
of the SET required in order to transfer an electron 
from the SET island into the lead or vice versa. Away 
from the anomaly region, the activation energy, around 
a sequential-tunneling peak, is given by 



A= \E{n,N) - E{n,N + 



(6) 



because the transport occurs due to the fluctuation N o 
TV + 1 at fixed n = (at Ng < N^) or n = 1 (at Ng > 
N^). According to Eq. ([5]), the activation energy changes 
linearly, on both sides of the peak, from the smallest 
value, A = 0, at the peak center, to the largest value, 
A = Ec, in the middle of each of the two adjacent CB 
valleys. However, as one approaches the anomaly region, 
the lower and the upper bounds of A change. 





FIG. 5. Activation energy A around a sequential-tunneling 
peak, (a) Upper panel: fragment of the stability diagram, see 
Fig. Ills, showing the working line (dash-dotted line) outside 
the anomaly region. Lower panel: the A vs. Ng dependence, 
corresponding to the upper panel, (b) Same as in (a), but 
with the working line in the anomaly region. The smallest 
value of A is (5. 



In this section, we consider the transport through the 
donor-SET system in greater detail and derive the results 



We illustrate how A is modified around the lower 
bound in Fig. [Sj The left panel shows the usual behavior 
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away from the anomaly region. The right panel shows 
the modification of A and the occurrence of a different 
lower bound, 5, in the anomaly region. Indeed, since the 
ground state is (1,A^) and the closest transport channel 
is (0, A^) o (0, + 1), then transport is possible only if 
the donor happens to be in an excited state, correspond- 
ing to charge n ~ 0, instead of the lowest energy state 
(n = 1). Thus, to the left of the peak center, we have 

A^E{0,N+1)-E{1,N), (7) 

whereas, to the right, we have 

E{0,N) ~ E{l,N). (8) 

This modification of the activation energy, as compared 
with Eq. ([6]), results in strongly asymmetric CB peaks at 
low temperatures. The peaks become sharper on the side 
which extends away from the center of the anomaly. The 
peaks maxima are still given by E{n, N) = E{n, + 1), 
as long as the temperature is sufficiently low (T <C 
UI2/ Ec)- The conductance at the peak maximum is sup- 
pressed as 

G oc exp(-(5/r), (9) 

where 5 is the lower bound of the activation energy A. 
From Eqs. ([7|) and (|8|), at the peak maximum, we obtain 

5 = [E{0, N) - E{1, A^)]U(o,M).i;(o,^+i) • (10) 

Equation (fTO|) is valid on the left-hand side of the 
anomaly, N- < Ng < {N+ + N-)/2. On the right- 
hand side, a similar derivation holds, with respect to the 
lower triple point, and S is given by 

6= [E{hN + l)-E{0,N+l)]\^^,j,^^E(,j,^,^, (11) 

for {N+ + N-)/2 < Ng < N+ . In the middle of the 
anomaly, 6 takes on the largest value, given by 

Here and throughout the paper we make the reasonable 
assumption that U12 < 2Ec- 

For simplicity, let us assume that Ng is measured with 
respect to the middle of the anomaly, i.e. N~ = —N^. 
Then, an explicit expression for S reads 

SiNg) = <5.nax (1 - ^) OiN^ - |iV,|), (13) 

where 9{x) = 1 for a; > and d{x) = for a; < 0. 
Equation (fT5|) is shown in Fig. HJ: as the solid line. 

In the special case — 0, the behavior of (5 is modified 
outside the anomaly region. Indeed, transport is possible, 
in this case, only due to hopping over the donor. Acti- 
vationless transport can occur only at the triple points, 
where two resonance conditions are met, 

r E{n,N) = E{n,N^\\ 

\ E{1,N) = E{0,N + 1), ^^^^ 



both for the upper (n = 0) and lower (n = 1) triple 
points. In this case, we find that Eq. ([T^ for 6 changes 
as follows 

S{Ng) = 5.nax (1 - ^) [2eiN+ - \Ng\) - 1] . (15) 
Equation (fTSj) is shown in Fig. Ht as the dashed line. 



2. Activation energy in the CB valley 

In Sec. nil A 1) we discussed the lower bound of the 
activation energy A, see Eq. ([S]). Here we focus on the 
upper bound of A and discuss the conductance in the 
CB valleys. We assume that the tunnel coupling is suffi- 
ciently small and, therefore, we discard cotunncling. 

We start with a position Ng, deep in the region of 
usual behaviour (for concreteness that with n ~ 0), where 
the activation energy is given by Eq. ^ and assumes 
A = Ec in the middle of a CB valley, say at Ng = 
N. It is important to remark that while the position 
Ng = N corresponds to the middle of a CB valley for 
the n — branch of resonances, this position is off the 
middle for the n — 1 branch. Thus, from the prospect 
of the n ~ 1 branch, the position Ng — N has a smaller 
activation energy because of the branch shift Ng ^ Ng + 
Ui2/2Ec- Therefore, if n would change from to 1, 
a transport sequence with a smaller activation energy 
would be possible. However, changing n from to 1 also 
costs an activation energy. The point at which the two 
factors balance each other is 

where we assume, as before, that Ng is measured with 

respect to the center of the anomaly. For |Afg| < N^ , the 
energy spent to change the donor charge is less than the 
gain in the activation energy of the transport sequence; 
and hence, A < Eq ■ The upper bound of A is given by 

6{Ng) =Ec-^(^l^ 0(Nf - \Ng\). (17) 

The positions of the minima shift linearly from Ng = N 
to Ng — N + U12/2EC as one traverses the upper-bound 
anomaly, Ng < Ng < N^ . 

3. Average occupation (n) 

In Eq. dSj, we expressed the linear conductance G{Ng) 
through the average occupation (n) . Here we give explicit 
expressions for (n) in different regimes. 

First, let us find the average occupation (n) of an iso- 
lated donor at thermodynamic equilibrium with a lead at 
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temperature T and chemical potential fi 
the general prescription, we have 

IPi + Opo, 

1, 

ed/T 



0. Following 





- Ipt + ] 


z = 






1 


Pt 


Pi = ^ 




1 


Po = 


Z' 



Thus, we obtain 



in) 



1 



1 



(18) 



(19) 



(l/2)exp(ed/T)- 

In the expressions above, po, p^, and pj^ denote, respec- 
tively, the probabilities for the donor to be in the empty, 
spin up, and spin down states; Z denotes the statistical 
sum. 

Next, we find the average occupation (n) of the donor 
in the presence of the SET, as used in Eq. (jS]). The 
general expression reads 



{n) = 



PnN 



npnN 

nN 
1 

: — Sr,e 

z 



E(n,N)/T 



(20) 



(21) 



where Z — J2nN Sne"^^"'^-*/^ is the statistical sum, s„ 
is the spin degeneracy on the donor {sq — 1 and si =2), 
and E{n,N) is given in Eq. ([T]). Equation (pO)) can be 
given a form similar to Eq. (|19p. 



(22) 



1 



l + (Z_/2Z+)exp[(erf-eO)/r]^ 



where e° 



N ' 



'-'12 

4Ec 
^-Ec(N 



is the middle position in the anomaly, 
Ng)^/T ^j^g statistical sum of the 
SET on the left side of the anomaly, and, correspond- 
ingly, Z+ = ^^e-^c(^-^»+^i2/2i5o)VT is the statisti- 
cal sum of the SET on the right side of the anomaly. We 
remark that, in the limit U12 <C T, Eq. ([^ reduces to 
Eq. (HH). 

Next, we consider the vicinity of a pair of triple points, 
as shown in the inset of Fig.|3}D. We assume T,Ui2 <^ Ec- 
The relevant charge states of the SET are N and -I- 1. 
Then, the factor Z^/Z+ in Eq. reduces to 



cosh [-^ 



cosh 



Eg 

T 



2Ec ) 



U12 , 



xe 



U12 



(23) 

Equation (|23|) together with Eq. (|22l) give the average 
occupation {n) in the vicinity of a pair of triple points. 

Around the upper triple point the charge configuration 
(1, iV -I- 1) can be omitted and the expression for (n) can 
be simplified further. 



in) 



1 



1 {l + e2Ec6NJT] 



,(<5ed-C/i2<5JVg)/T 



(24) 



where SNg = Ng — N — 1/2 and 5ed = — U12/2 are, re- 
spectively, the dimensionless gate voltage and the donor 
energy measured from the upper triple point. Similarly, 
around the lower triple point the average occupation (n) 
can be written as follows 



(n) 



1 



1 



^2Ec&NJT 



_ _ (25) 

where 5N g = Ng- N -1/2- Ui2/2Ec and Se^L = ^d + 
Ui2/2 — Ul2/'^Ec are, respectively, the dimensionless gate 
voltage and the donor energy measured from the lower 
triple point. 

Note that at the upper triple point the average charge 
is (n) = 1/2, whereas at the lower triple point it is (rt) = 
4/5. This asymmetry is due to the spin degeneracy on 
the donor. 



B. Kinetic theory of conduction 



In Sec. nil Al we discussed the activational conductance 
through the donor-SET system and related the activation 
energy A to the electrostatic model of Eq. ([T]). Such an 
analysis helps to understand the qualitative behaviour of 
the conductance G{Vg) in the anomaly region. A more 
rigorous consideration of G{Vg), would involve inclusion 
of a prefactor in front of the exponential dependence in 
Eq. (O. For large A/T (i.e. for small temperatures), 
such a prefactor is of little relevance, since most features 
in Gi/Yg) come from the dependence of A on Vg, and 
are enhanced by the fact that A appears in the expo- 
nent. At higher temperatures, the dependence of the 
prefactor on Vg can be comparable to the dependence 
of exp(— A/T) on Vg. It turns out (see below) that the 
analysis made in Sec. IIII Al is valid at T ^ Ui2/Ec'- In 
the rest of the paper, we give a more rigorous treatment 
of the problem. We complete the electrostatic model of 
Eq. ([T]) with kinetic-energy and tunneling terms, and give 
a rate-equations description of the transport problem. 
We analyze different temperature regimes and find that, 
at r ^ U12, the rate-equations description is identical to 
a description in which the actual (interacting) system is 
replaced by a circuit of resistors. We employ such a re- 
placement as an approximation at any temperature and 
illustrate its accuracy for the parameters relevant to the 
experiment. 



1. Generalities on the Hamiltonian and sequential 
tunneling theory 

We assume that transport through the SET is due 
to sequential tunneling , Si i^s , 34 qj-^^qj. calculate the 
linear conductance, we rewrite the model in Eq. ([T]) in 
second-quantized form and complement it with tunnel- 
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ing terms. The donor Hamiltonian takes the form 

Hd^ <id ^ rf^-dcr + U^n^n^, (26) 
<T=t4 

where d\ {da) is the creation (annihilation) operator of 
an electron of spin a = {'f , 4-} on the donor. The donor 
occupation is given hy n — + n^, where = d^^d^ is 
the occupation of one spin species. The on-site repulsion 
energy Uoa on the donor is assumed to be large enough in 
order to exclude the occupation n — 2. The SET island 
is described by 

Hd = J2 ^kflfka + Ec{N- Ngf , (27) 

k{7 

where Sk is the energy spectrum in the island, fl^ (fka) 
is the creation (annihilation) operator of an electron in 
state {k,a) on the island, and N = J2ka flafko- the 
operator of the number of electrons on the island. The 
single-particle energy spacing in the island is assumed 
to be vanishingly small (much smaller than the temper- 
ature). The electrostatic and tunnel coupling between 
donor and SET island is given by 

HdD - Uu (N - Ng) *i2 J2 ("^yka + h.c.) , 

(28) 

where ti2 is the tunneling amplitude between donor and 
SET island. The leads, i.e. the source (L) and the drain 
(R), are described by 

-ffleads = 'l2^p4pa^lP'^' (29) 

Z = L,fl per 

where is the dispersion relation in the leads and cj^^^ 
(c(po-) is the creation (annihilation) operator of an elec- 
tron with momentum p and spin a in lead I — L,R. The 
tunneling between the donor and the left lead is described 



by 

HdL = tL Y^idULpa + h.c), (30) 

pa 

where is the corresponding tunneling amplitude. The 
tunneling between the SET island the leads is described 

by 

Hdlr Y^fLcipa + h.c), (31) 

I kpa 

where Vi is the tunneling amplitude between SET island 
and lead I = L, R. 

Next we calculate the sequential-tunneling rates for 
each of the tunneling terms in Eqs. (|28|) . (pO| . and (|3T1) . 
The sequential tunneling between donor and SET is de- 
scribed by the rates 

yyiN,oN+i — 

X fied + Ui2{N - Ng) - Ec {2{N - Ng) + 1)) , 

yyQN+i,iN — [t 

- f{ed + U^2[N ~ Ng) ~ Ec {2{N - Ng) + 1))] , 

(32) 

where Tr = irvjj |ti2|^, with being the density of 
states in the SET. The different prefactors in the direct 
and reverse rates come from the spin degeneracy on the 
donor. The sequential tunneling between donor and lead 
is described by the rates 

WiN,0N = -f^f {ed + Ui2{N - Ng) - AfiL) , 

WoN,iN = ^ [1 - / (q + Ui2{N - Ng) - A^il)] , 

(33) 

where Tl = t^v \ ti}^ with v being the density of states in 
the lead, and A/j,; = ^i — ^, for I = L,R. The sequential- 
tunneling rates of the SET at a fixed n read 



WnN+l,nN - ^ W^j,+ ,^oN ^ Y (^^' " '^^^ {N - Ng + 1/2) - Ui2n) ., 



l=L,R I 



Wr,N,nN+i= V W'^^ +^ ^ Y (Ui2n + 2Ec {N ~ N g + I /2) - Ain) , (34) 

l=L,R I 



where gi — An'^iyvi:, \Vif is the dimensionless conductance (in units of e'^/irh) of lead I and Q{E) — 
E/ [l-cxp{-E/T)]. 



With the help of the rates in Eqs. (|32p . ([33]). and (p4|) . we are able to calculate the linear conductance of the 
donor-SET system in different regimes. 
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2. Two triple points (Ui2,T <^Ec) 

We consider the vicinity of two triple points, as shown in the inset of Fig. |4}3, and write down the rate equations 
for four charge configurations, (OA^, OiV + 1, lA^, lA^ + 1), 

(WoAr+l,OAr + W^lW.Ow) PON — Wqn^oN+iPoN+I + Won,1NPiN, 
{WoN,ON+l + WiN,ON+l + M^lAT+l.OW+l) PON+1 = Won+1,ONPoN + M^OAf+l.lA^Pl^ + Won+i,in+iPin+i, 
{Wqn,1N + Wqn+i^in + M^lAr+l,lAr) PiN — WinmnPoN + Win,0N+iP0N+1 + M^lAf,17V+lPl7V+l , 

(WoAT+l.liV+l + WiN,lN+l) PlN+1 — WiN+l,ON+lPON+l + Win+1,1NPiN- (35) 



These equations are comphmented by the normahzation 
condition poN + PoN+i + Pin + PiN+i = 1- The cur- 
rent can be evaluated, for example, at the right tunnel 
junction. 



I = - \e\ (Wg^ o^_|_iPo7V+l + W^/Ar,lAr+lPlAf+l 



-^ON+lfiNPON - W^lAr+iaTvPlAr) 



(36) 



\e\dI/dA^i 



The linear conductance is obtained as G = 
in the limit — > 0, where A/i — ~ PR- 

We focus in the following on the case = and 
assume, for simplicity, that F^, ^0. This is the case 
of interest with regard to Eq. ([3]). Using the rates in 
Eq. (|34p , we calculate the conductance through the SET 
for a fixed charge on the donor and recover the well- 
known expression 



9LgB. 



h Ql + gRS\nh{x - xnY 



(37) 



where x — 2EcNg/T . The position of the peak is given 
by XN = 2Ec{N + 1/2) /T for the donor in state n = 
and by a; AT = [2Ec{N + 1/2) + U12] /T for the donor in 
state n = \. In the limit F^ — !> 0, the donor changes state 
seldom and the current as a function of time represents a 
telegraph noise. Such telegraph noise is observed in our 
samplesSS, for a range of the gate voltage, consistent with 
the picture given here. 

The probability for the donor to be in state n = is 
given by po = 1 ~ and the probability to be in state 
n = 1 is pi = {n). Here, {n) is the average occupation 
on the donor, considered in Sec. IIII A3I The linear con- 
ductance is the average over the two realizations of the 
donor occupation. 



e 9L9R 



X — XnN 



h gL + 9r ^ smh(a; - XnN) 



(38) 



where XnN = \2Ec{N + 1/2) + J/i2n] /T. We note that 
the only way how correlations between the donor and 
the SET charges enter in Eq. ([55]) is through the av- 
erage occupation (n), which contains information about 
the feedback of the SET state on the donor occupation 
probability. For T ^ Ec^ we sum the right-hand side of 
Eq. dSH]) over N and obtain Eq. ©. 

The periodicity of CB oscillations is broken in the 
anomaly region. The CB peaks occur at positions f^g = 



-I- 1/2 to the left of the anomaly and at positions 
TVg = iV -h 1/2 U12/2EC to the right of the anomaly. 
The question arises: How is the shift of CB peaks dis- 
tributed in the anomaly region? Let the CB peaks in the 
anomaly region occur at position 



Ng = N+^ + ^{Ng), 



(39) 



where the phase shift f{Ng) changes from to Ui2/2Ec 
as one crosses over the anomaly. We find that there exist 
three temperature regimes in which ip{Ng) has different 
behaviour. 

At low temperatures, T ^ U12/EC, as we already 
mentioned in Sec. I A, the CB peaks in the anomaly 
region have asymmetric shape. The conductance decays 
exponentially on both sides of the peak. However, the ex- 
ponents of this decay are not the same on both sides, see 
Fig. m The peak is sharper on the side that faces away 
from the center of the anomaly. In the low-temperature 
limit (T — !■ 0), the position of the peak is determined by 
the smallest activation energy. The dependence of on 
Ng is a step function. 



2Ec 



(40) 



with the jump occurring in the middle of the anomaly. 

At intermediate temperatures, U12/EC <C T <C C/12, 
the activation energy A, entering in Eq. ([5]), can be con- 
sidered to be constant in a region of Ng where A varies 
weakly, see Fig. [S] The shape of the CB peak is deter- 
mined by the prefactor of the exponential dependence 
in Eq. ([5]) and resembles the sharp angle of a right tri- 
angle: the conductance increases linearly, reaches max- 
imum, and then drops abruptly. The sharp side of the 
peak faces the center of the anomaly. When we speak of 
the peak position, we understand the position given by 
the first moment of Ng, 



_ jNgG{Ng)dNg 



jGiNg)dNg 



(41) 



This position differs from A^ -I- 1/2 by the shift 



12 



6Ec 



(W^-IN,!)^ (42) 
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Here we neglected small corrections proportional to tem- 



perature in the vicinity of Ng — 



and a broaden- 



ing of the step function over a small region of Ng, given 
by l^^gl < {T/Ui2)N+. We remark that, if, instead of 
Eq. (I4ip . we use the position of the peak maximum in 
order to talk about the shift (/s, then Eq. (|42|) is modified 
as follows, 



f ■ 



Ul2 



^ + ^]HK-\Ns\)- (43) 



At high temperatures, T 3> U12, the dependence of (n) 
on Ng is weak. Substituting Eq. ([55]) into Eq. (gT]) and 
taking (n) out of the sign of integration, we obtain 



12 



2Ec 



[n). 



(44) 



AN, 



S/Ec 



T/Ec 



T/Ui 



^ {Ui2/Ec)5 5 T 

FIG. 6. The width of the CB peak in the anomaly region as 
a function of temperature. 



An expression for (n) in this regime is given in Eq. (|19p . 
We remark that this result is consistent with the mean- 
field approach, in which one takes the average of n in the 
system Hamiltonian and then calculates the SET con- 
ductance. The resulting effect is a shift of the conduc- 
tance peaks by an amount given in Eq. (|44p . The shift 
of CB peaks is distributed over a region of Ng which ex- 
tends outside the anomaly. In practice, such a shift can 
be detected only by comparing the shifts of neighbour- 
ing peaks. For this purpose, the following expression is 
useful. 



3. One triple point (T < U12 , Ec ) 

We now focus on the upper triple point, retaining three 
charge states, (OiV, 0A^+ 1, IN), see Fig. [5] Around a sin- 
gle triple point, Eqs. (|55|) and (pS)) allow a more explicit 
consideration. Quite remarkably, the end result of this 
lengthy calculation is an expression which can also be ob- 
tained by assigning an effective resistance to each tunnel 
junction and replacing the studied system by a circuit of 
resistors, which is in the spirit of Ref. [3^ 



dip 



aUr. 
2T 



(45) 



Furthermore, the sample may also contain donors which 
are situated further away from the SET, and therefore, 
also from the top gate. For such donors, the quantity 
in Eq. (H5|) can be neglected, since both a and U12 are 
small. The leading-order effect comes from Eq. (|44)) and 
represents a shift ipQ which is constant as a function of 
Ng. We expect (po to have activational dependence on 
temperature, tpo oc exp (— Ajmp/r), with Ajmp being the 
activation energy of the distant donors. 

Next, we discuss the width of the CB peak as a func- 
tion of temperature. Outside the anomaly region the 
width is proportional to the temperature, AiVg T / Ec- 
At high temperatures, T ^ U12, this results persists in 
the anomaly region. At lower temperatures, however, 
the width of the CB peak is considerably larger in the 
anomaly region than it is outside. We find that the 
width of the CB peak, depends on the position Ng in 
the anomaly, Ng' < Ng < N^. This dependence en- 
ters through the dependence of 6 on Ng, see Sec. I A. We 
summarize our results in Fig. [5] At intermediate temper- 
atures, {Ui2/Ec)5 ^ T ^ (5, the width is temperature- 
independent and given by ANg ~ 5/ Ec- At low temper- 
atures, T -C {Ui2/Ec)S, the width is again proportional 
to temperature, but with a different proportionality co- 
efficient, ANg - T/U12- 



Source 
— 



Donor 



rR 



SET 
island 



Rr 



Rl 



Drain 
I — 



FIG. 7. Circuit of resistors, representing the donor-SET sys- 
tem at low temperatures. 



The donor-SET system can, thus, be represented by 
the circuit in Fig. [71 The resistances associated with 
each of the tunnelling terms (see rates in Eqs. (|32|) . ([33|), 
and ([Ml) ) are given by 

1 

tr T T 
1 

TL T T 

1 _ e2 ; _ e2 ; 

(46) 

where all the quantities are evaluated at thermodynamic 
equilibrium. We note that each of the expressions in 
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Eq. p6)) has the form of the Einstein relation, 



G = 



T 



(47) 



where the role of the diffusion constant D is played by 
the product W p with appropriate indices. The total con- 
ductance follows from adding the resistances in Fig. [71 



i ^ „ RL{rL+rR) 
G RL+rL+rR- 



rn 



To give explicit expressions for the resistances in Eq. ([i^l , 
we measure energy with respect to the charge configura- 
tion (0,-/V) and obtain 



PON — 

p^^ ^ 1 2e-(^^--^-^^»)/^, (49) 
where Z is the statistical sum given by 

Z = 1 + 2e-(^^<^-^i^*^«)/^ + e^^^^^'ol^. (50) 

For the notations of bt and 5Ngi see expressions in the 
text below Eq. (|24p . The final expressions for resistances 
read 



TL 



Ri 



h ZT 

^4f^ 
h ZT 
^4f7 



Ui2SNg)/T 



e2 gi 2EcSNg ■ 



(51) 



In the rest of this section, we derive Eq. ([4]). We set 
gL = and assume gR ^ ^i/T. In this case, the linear 
conductance is determined by the resistance of the donor 
tunnel junctions, 1/G = tl + vr, see Eq. (pS]). The linear 
conductance is largest at the triple points (or close to the 
triple points within a distance oc T) . The CB oscillations 
show, therefore, a resonance in the envelope function of 
the peaks when the working line of the device intersects 
either the upper or the lower row of triple points. The 
questions of interest, here, are: What is the shape of the 
envelope function and how does the maximal value of 
conductance at the resonance scale with temperature? 

The envelope function can be determined by, first, find- 
ing the positions of the maxima of the CB peaks, and 
then, evaluating the conductance at these points. With 



the simplifying assumptions we made above, the posi- 
tions of the maxima are determined by solving 



d 



dSN, 



{tl + 7-_r) = 



(52) 



with respect to 6Ng. This task, however, makes the end 
result cumbersome and difficult to discuss. We aim, in- 
stead, at finding a representative curve for the envelope 
function, i.e. a curve which is correct by order of magni- 
tude. We find that, for Fj^ < F^, the solution of Eq. ^ 
can be approximated as 5Ng w 0. Setting 5Ng = in the 
expressions for and tr in Eq. ([51]) and using Eq. ^ 
we obtain 



1 hT T, 



G 2e2 TlTi 



l-t-e 



IT . 



1 



(53) 



which derives Eq. provided we also relate to Ng 
using the "working line" relation. 

Focusing on the lower triple point, i.e. retaining only 
the charge states (OiV -|- 1, \N ^ IN + 1), we derive sim- 
ilar expressions to those above. In particular, Eq. ([55]) 
becomes 



1 hT T, 



(4 -t- e^'"/^) (e-*^''/^ + 1) , (54) 



G 4e2 TlTi 



where Scd is defined in the text below Eq. ([25]) . 

We conclude this section by remarking that (i) the en- 
velope function of the CB oscillations resembles the reso- 
nance transmission of the donor and (ii) the temperature 
dependence of the maximum of the envelope function is 
Gocl/T. 



4- Discussion of a "sum- of -resistance" formula for 
arbitrary case 

Here, we derive an approximate expression for the lin- 
ear conductance by extending the circuit approach of 
Sec. nil B 31 to arbitrary temperatures. We remark that 
the circuit approach is a quick and computationally- 
inexpensive way to obtain an insight into the behaviour of 
the linear conductance. The approximation made when 
replacing the system under consideration by a circuit of 
resistors consists in neglecting non-local correlations, i.e. 
correlations extending over several tunnel junctions. We 
compare the circuit approach against a numerically exact 
evaluation of the conductance at the end of this section. 
Below, we derive expressions for the resistances of the 
resistors in Fig. [7| 

We consider a weak deviation from equilibrium, for 
which it is not necessary to correct the distribution func- 
tion. The current flowing from the lead I = L, R into 
the SET island originates from the fluctuation (n, N) o 
{n,N 1) and can be written as 
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^' = ~^E / ds'{pnNfie)[l-fis')]-pnN+l[l-fmi^')} 



X(5(e - e' + A^; + E{n, N) - E{n, N + 1)), 



(55) 



where pnN is the equiUbrium distribution function, given 
in Eq. (PT|) . and A/i; 0. Then, the resistance Ri asso- 
ciated with junction I is found from Eq. ([55]) to be 



J2pnNQ{E{n,N)-E{n,N + l)). (56) 



In the same fashion, we write down the current flowing 
from the SET island to the donor, considering the fluc- 
tuation (0, N +1) ^ (1, N), 

iR = ^ [ '^^ {pin[1 - /(e)] - 2pow+i/(e)} 

N •' 

xS{e + ApD+E{0,N + l)-E{l,N)). (57) 

Here, the factor of 2 in front of poN+i originates from the 
spin degeneracy on the donor and we assumed that the 
SET island is at an elevated chemical potential Afijj — )■ 0. 
The resistance rn associated with this junction is, then, 
found to be 



1 4r 



rR h T 



N 



(58) 

Similarly, the resistance associated with the junction 
between the donor and the left lead is obtained from 



rL fi T 



Y,P0Nf{E{l,N)~E{0,N)). (59) 



N 



With the help of Eqs. and (HH), the conduc- 

tance G is given by Eq. (jlSj). 

In Figs.[3t and[3ji, we reproduce the measurement data 
of Figs. andHjD, respectively. The parameters we have 
used in Eqs. (gll), (HH), ([58]), and §9^ are listed in TableHl 

Next, we compare the circuit approach against an ex- 
act numerical evaluation of the conductance. In Fig. |S1 
we show Figs. ^ and[3Jl (solid line). For the same set of 
parameters, we evaluate the conductance exactly (dots) 
by solving the Pauli master equation numerically. In or- 
der to do this, we generalize the set of equations ([55)) 
to include all charge configurations {n,N). Then, we re- 
tain a sufficiently large number of charge configurations, 
around N ?» Ng, to guarantee convergence, especially for 
T > Ec- We find good agreement between the two ap- 
proaches for the parameters relevant to the experiment, 
see also caption of Fig. [S] 



IV. CONCLUSIONS 

In conclusion, we have used a SET as a highly sensi- 
tive device to study single dopant atoms in Si. A simple 



theoretical model allowed us to explain the anomalous 
behaviour of the linear conductance seen in the experi- 
ment and to extract values of the coupling strength for 
the dopant atom and the SET. Our results can be used 
to assess the feasibility of using a SET as a means to 
manipulate and read out single dopant atoms in Si. 
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Appendix A: Stability diagram 

The stability diagram, see Fig. |3l3, shows regions in 
the two-parameter space {Ng,ed), where the lowest-in- 
energy charge configurations of the donor-SET system 
is definite. The solid lines separate the regions of dif- 
ferent charge configurations, and hence, represent lines 
at which the charge configuration is degenerate. Using 
the electrostatic energy in Eq. ([1]), we describe here the 
stability diagram quantitatively. 

The upper vertical lines in Fig.HjD correspond to degen- 
eracies of the type E{0, N) = E{0, N +1). The position 
of these degeneracies is given by 



(Al) 



where N = 0, ±1,±2... takes on integer values. The 
inclined lines with positive slope correspond to E{0, N) = 
E{1, N) and are given by 



ed^Ui2{Ng~N). 



(A2) 



The inclined lines with negative slope correspond to 
E{0, N + l)= E{1, N) and are given by 

ed^Ec-{2Ec-U,2)iNg-N). (A3) 

The lower vertical lines correspond to E{1,N) ~ 
E{1, N + 1) and are given by 



1 , Uu 



2 2Ec 



(A4) 
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TABLE I. Parameters used in reproducing the experiment in Fig. [3] Samples 1 and 2 denote, respectively, the samples of 
Figs. [3^ and[3l3. The parameters in each column are used in Eqs. (|48p . (|56p . (|58|l . and (|59|l to reproduce the measurement data 
at the corresponding temperature. The result for Sample 1 (2) is shown in Fig. |3j; (Fig. |3ji) . Parameters a and jS determine the 
position and orientation of the working line with the help of the relation = Ui2/iEc — aEc{Ng — /3). The relation between 
Vg and Ng is conveniently parametrized as follows, Vg = Vg{0) + {e/Cg)Ng, where the aim of the offset Vg{0) is to restrict Ng 
(and hence A'^) to small numbers. 




1.15 1.2 1.25 1.3 1.35 1.45 1.5 1.55 1.6 



Vg (V) Vg (V) 

FIG. 8. Comparison of the circuit approach against an exact numerical evaluation of the conductance. The solid lines represent 
the result of the circuit approach, whereas the dots show the conductance evaluated by solving the Pauli master equation 
numerically, (a) Conductance calculated for the same parameters as in Fig.[3t. The two approaches coincide identically in this 
case, because the circuit approach is exact for the donor-SET model in the limit Fj — i- 0, (/ = L, R). In particular, Eq. (|38p can 
be obtained from Eqs. (I48|) . (|56p . (|58p . and (|59p . (b) Conductance calculated for the same parameters as in Fig. [3ji. For the 
low temperature trace (T = 0.9 K), the circuit approach deviates from the numerically exact result by a relative error e < 36%. 
However, the deviation takes place far in the tails of the peaks where the conductance is small; hence no visible difference on 
the scale of the plot. For the high temperature trace (T — 4.2 K), the relative error is £ < 7% and the deviation is visible in 
several CB valleys. 



The upper triple points correspond to the double con- 
dition E{0,N) = E{1,N) = E{0,N + 1) and have the 
coordinates 



Ed = 



12 



(A5) 



Similarly, the lower triple points are found from E{1, N - 
1) = E{1, N) = E{0, N + 1) and have the coordinates 



Ul2 

2Ec' 



fid 



r/2 



Ul2 

2 ' 2Ec' 
5l and (IA6I) one finds that the 



(A6) 



From comparing Eqs 
middle position in the anomaly (middle between upper 
and lower triple points on e^-axis) is given by 



AEc 



(A7) 



Similarly, the width of the anomaly (distance between 
upper and lower triple points on e^-axis) is given by 



^12 

2Ec' 



(A8) 
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Using the equation for the working hne of the device, 
Cd — —aEcNg + const, we relate A^^, see Fig. to 
the parameters of the stabihty diagram. Considering the 
edges of the anomaly, we write 

e° T ^ = -aEcNf + const. (A9) 

Together with Eqs. (jA7[) and (jA8l) . the latter equation al- 
lows one to express N'^ in terms of the parameters C/12, 
Ec, a, and the offset (const). It is convenient for a prac- 
tical purpose, to parametrize the offset by a parameter 
/3, such that the equation of the working line is written 
as follows 



^d = -^-aEc{Ng-l3). (AlO) 



The values of /3 can be restricted to the interval < /3 < 1 
by introducing an offset to the gate voltage Vg, 

V,^VgiO) + '-^, (All) 

where Vg{0) is the gate voltage offset. With this choice, 
the charge configuration with = lies close to the 
center of the anomaly, and the parameter /3 gives the 
position of the anomaly center, Ng = (3. The anomaly 



edges are, then, given by 

".*-^-2Sfe('-^)' <*'^) 

In the limit a ^ 1, when discussing envelopes of CB os- 
cillations, one may dispense with the term /3 in Eq. (jA12|) . 

Finally, we remark that, for a fixed value of a, the 
width of the anomaly is largest at U12 = Ec, see 
Eq. (|A12[) . and that the anomaly region shrinks to a point 
when U12 approaches its largest value, U12 — 2Ec- The 
absence of the gapped region for U12 — 2Ec is related 
to the fact that the electrostatic model in Eq. ([1]) can be 
rewritten as follows 

E{n,N) = (ed-Ec)n + Ec{n + N -Ngf, (A13) 

where we used the identity — n. The first term on 
the right-hand side in Eq. (jA13l) can be interpreted as a 
kinetic energy term, and thus, can be included into the 
kinetic-energy term of Eq. ([27l) as an additional energy 
level in the SET island. The second term represents the 
usual Coulomb energy of an island, with the difference 
that the donor appears as part of that island too (the to- 
tal charge being n+N). It is important to note, however, 
that this analogy to the usual SET (and hence a seem- 
ing disappearance of the donor from the problem) goes 
as far as the electrostatic energy is concerned and spin is 
ignored. The donor may still leave visible traces in the 
transport, even sX U12 = '2Ec, provided both and Tji 
are nonzero. As a matter-of-fact, the low-temperature 
trace of Fig. [2^ suggests precisely this scenario. 
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